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•^ ■ Abstract 

r-| ■ The standard description of material media in electromagnetism is based 

^ Oh | on multipoles. It is well known that these moments depend on the point of 

reference chosen, except for the lowest order. It is shown that this "origin 
dependence" is not unphysical as has been claimed in the literature but forms 
only part of the effect of moving the point of reference. When also the com- 

ff^ I plementary part is taken into account then different points of reference lead 

^ ■ to different but equivalent descriptions of the same physical reality. This 

Q , is shown at the microscopic as well as at the macroscopic level. A similar 

VO I interpretation is valid regarding the "origin dependence" of the reflection co- 

^ ■ efficients for reflection on a semi infinite medium. We show that the "trans- 

^ . formation theory" which has been proposed to remedy this situation (and 

• 1— I I which is thus not needed) is unphysical since the transformation considered 

J>L,' does not leave the boundary conditions invariant. 

^: 
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1 Introduction 

In classical electrodynamics! 1 1 material media are modelled at the microscopic 
level as an ensemble of stable building blocks (atoms, ions, molecules, . . . ) with 
certain microscopic charge (r/) and current (j) densities. The microscopic fields 



(e,b) then obey the microscopic Maxwell equations: 

db 

Vxe = -— V-b = 

^ b dene - ^ 
Vx — = — ^ + ; V-eoe = 7] 

Mo ot 

With a purely classical model the microscopic sources may be written as: 

r\=Y^qkl5{r-Rki) j = Y,1klRklSir-Rki) (2) 

where the index k labels the building blocks and I the point charges (with po- 
sition Rj^i and velocity i?^/) within a building block. The macroscopic Maxwell 
equations then follow by a suitable averaging procedure. Due to the linearity 
these equations are similar to the ones in (HI) with (e, b) replaced by the macro- 
scopic fields {E,B) and with the microscopic sources replaced by their averages 
((77), {])). The microscopic source densities Q can be rearranged and expressed 
in terms of the (microscopic) electric and magnetic multipole moment densities. 
Usually only the 2 lowest order multipole moments are taking into account. After 
averaging, these lead to the usual macroscopic source densities: 

- - dP - (3) 

{j) = j+ —+yxM 

where (p,7) are the unbound charge and current densities and P and M the 
polarization and magnetization densities. In order to avoid the complexities asso- 
ciated with moving media, in what follows we will consider only non-moving me- 
dia. Substituting Q in the averaged ([T]) and defining D = £qE-\-P and H = -^-—M 
one finds the standard macroscopic Maxwell equations: 



- dB 
dt 

VxH = ^+J VD = p 

dt 



(4) 



If higher order multipole contributions are taken into account then additional 
terms show up in ©. By adapting the definition of iD,H) accordingly equations 
^ remain valid. However it should be noted that at a boundary between 2 media 
the jump conditions which follow from Q have to be adapted in order to take into 
account these additional multipole contributions. 



It is well known that the multipole moments are not unique and more precisely 
they depend on the choice of the reference point (e.g. the centre of mass) used for 
calculating the moments, except for the lowest order non vanishing moment. E.g. 
the total charge of an ion q^ = E/ ^« is independent of the reference point but 
the electric dipole moment is not, except if ^^ = 0. For insulating materials the 
dipole terms are unique and since higher moments are usually not considered the 
non-uniqueness of the higher moments is often of no consequence. 

In Q material media are treated up to the electric octopole and the mag- 
netic quadrupole terms. The emphasis is on optical phenomena (in insulating 
media) and in particular Fresnel reflection coefficients are calculated for a semi- 
infinite medium taking into account these quadrupole and octopole effects. There- 
fore the authors are faced with the problem of the non-uniqueness of these higher 
order moments. As a result of lengthy calculations they find that in general those 
reflection coefficients also depend on the chosen reference point (it is called "ori- 
gin dependent", but this should not be mistaken for the origin of the laboratory 
frame, but as a reference point within an atom.) These results are then discarded 
as being "unphysical" and using a "transformation theory" they adapt the "stan- 
dard theory" so that "origin-independent" reflection coefficients are obtained. Al- 
though at first sight the reasoning may seem tempting no indication is given at 
what point exactly the "standard theory" became "unphysical". After some con- 
sideration we became convinced that there is actually no problem with the "stan- 
dard theory" and that on the contrary the "origin dependence" as calculated in ^ 
has no physical meaning since only part of the effect of changing the position of 
the "origin" has been taken into account. In fact when choosing another reference 
point not only the higher order multipole moments do change but so does the po- 
sition of these elementary sources. If both effects are taken into account then the 
microscopic charge density remains unchanged as expected. This is explained in 
more detail in section 2. In sections 3 and 4 the same idea is worked out in the 
macroscopic domain and in section 5 the effect on reflection coefficients is con- 
sidered. It will become clear that there is no need for a "transformation theory". 
On the contrary, since the "standard theory" and the "transformation theory" give 
different results, at least one of them must be unphysical. In the last section 6 
it will be shown that the "transformation theory" is unphysical since it does not 
leave the boundary conditions invariant. 



2 Microscopic charge and current densities 

In order to simplify the notation we consider only a single building block. Then 
up to quadrupole order the microscopic charge density is written as: 

T] = qdir-R) - p,d,d{r~R) +q,jd,j5ir-R) + ■ ■■ (5) 

where the electric multipole moments are defined aslHI: 

^^"^ = ^L^rr^ (6) 

n. ^ 

In ^ we use {q.pi.qij) as a shorthand notation for the 3 lowest order mul- 
tipole moments, dt for the derivative with respect to xt, the cartesian components 
of F, R for the position vector of the reference point and 7/ for the position vector 
of charge I with respect to this reference point (thus Ri = R + Ti, where compared 
with Q the index k has been omitted). If we choose another reference point 
within the building block via R = R +d and F/ = r'i — d then in general different 
multipole moments are obtained: 

Pi = Y, ^i^ii = II <?/ i^'ii -di)= Pi - qdi (^Y^ 

and similarly for the quadrupole moment: 

(lij = 4ij - 2 i^iPj + Pi^j) + 2^didj (8) 

On the other hand these modified moments must now be placed in the new 
reference point. With a Taylor series we write: 

d{r-R)=d{r-t -d) 
_ _ 1 _ (9) 

= 5 (F - i?') - didi8{r - r) + -didjdij5{r -R') + --- 

Inserting Q (EJ) and ^ into ^ and retaining only terms up to quadrupole 
order one finds a similar expression but referred to the new reference point: 

r] = qdir-R') -p'^d,d{r-R')+q'ijdijd{f-R') + ■■■ (10) 

Therefore the microscopic charge density (r/) is as expected independent of 
the choice of reference point within the building block, although the multipole 
moments themselves do change according to © and ©. The interpretation is 
straightforward: suppose q ^0 then moving the reference point over —d changes 



the dipole moment with qd but the monopole q must be moved with the reference 
point and this constitutes an additional dipole moment exactly cancelling the first 
change. 

A similar conclusion can be drawn with respect to the microscopic current 
density. We define the magnetic multipole moments by: 

v^"^-(;^i:^/^r^(^/xrz) (11) 

with in particular for the magnetic dipole moment {n= 1): 

ini = -l^qieijkXjiXki (12) 

where etjk is the Levi-Civita tensor. This dipole moment transforms according 
to: 

mi = m'i - -eijkdjp'y- (13) 

Up to electric quadrupole and magnetic dipole order, the microscopic current 
density is then given by: 

ji = Pid{f-R) -qjidjd{r-R)+e,jkmkdjd{f-R) (14) 

Substituting dTJ) dH) (fT3t and Q into (fT4l) again results in a similar expression 
for the microscopic current density with respect to the new reference point: 

ji = Pid{r-R') - qjidjd{r -R') +eijkmkdj5{r- r') (15) 

We consider now as an example the reflection of a plane wave impinging on 
a semi infinite medium as was treated in [2J . In principle the problem could be 
solved at the microscopic level and this would result in definite values for the 
charge and current densities in each atom of the medium and thus for the global rj 
and j. We could then describe these resulting charge and current densities using 
multipole moments up to some chosen order and as shown this representation is 
independent of the choice of reference point, although the moments themselves 
are "origin dependent". This shows that already at the microscopic level the "ori- 
gin dependence" of the multipole moments has no physical consequence supposed 
they are treated in a consequent way, meaning that the multipoles should be placed 
at their particular reference point. 



3 Macroscopic charge and current densities 

Since the macroscopic charge and current densities ({rj), {j)) are obtained by av- 
eraging the microscopic densities, this conclusion must extend to the macroscopic 
level. The proper way of averaging is over an ensemble [5| with a probability 
density f{Rki,Rki\t) which depends on the positions of all charges and their ve- 
locities. Alternatively the probability can be expressed as a function of the atomic 

positions (Rk,Rk) and the internal coordinates irki.Tki). Since the probability of 
finding the system in a particular configuration is independent of the choice of the 
atomic reference point it follows that: 

f(Rk,%,rklJkht)=f'(R'kXj'kl,rM;t) (16) 

where we have already taken into account that the Jacobian of the transforma- 
tion (Rk = R k + d,Tki = f'ki — d) equals unity. If we take the sum over all building 
blocks of the first contribution in © and take the average over this distribution 
function, we find the unbound charge density p in (jSj): 

p = (^qk8{r-Rk))=Y,Qk / 5{r-Rk)fk{Rk;t)dRk = Y,^kfk{r,t) (17) 

k k ■^ k 

where /^(r; t) is the probability density for finding building block k at F. Again 
for simplicity we consider now only the contribution of one particular kind of ions 
and (llTt can then be written as: 

p=qh{r-t) (18) 

where /i (r; t) is the probability density for finding 1 ion of this kind at F. 
Note that (i) since we excluded moving media -^ = and (ii) p should be sup- 
plemented with a contribution from the free electrons; but since these have no 
internal structure they are not relevant for our subject. 

If we choose another reference point (shift of the origin) then in general an- 
other unbound charge density (due to the ions) is found: 

p' = qf\{r,t) (19) 

Due to (IT6b the probability densities are related by fi{r,t) = f\{{j — d)\t) 
and therefore p (F) = p' {? — d) . After developing we find: 

p = p'- d,d,p' + ^-didjdijp' + • • • (20) 

This result is completely similar with that found for the corresponding mi- 
croscopic quantities where the 5-functions in the latter are here replaced by the 



probability density function. The same reasoning can be followed for the higher 
order terms. For the polarization we find: 

P={Y,Pk^('^-^k))= P{n)f2{rji;t)dri (21) 

k '' 

where the integration is over all internal coordinates and a similar expres- 
sion holds for the shifted origin where the distribution functions are related by 
fi^Xi^ht) = f'2^~dyi'^^)- Taking into account Q and ( IT9b we then find: 

P,{r)=P',{r-d)-dip'{r-d) (22) 

where we have also taken into account that the integration of f'2 over the 
internal coordinates yields f\ and where we have omitted the time dependence. 
After developing we find: 

Pi = P'i - dip' + didjdjp' - djdjP'i + ■■■ (23) 

For the quadrupole density we find in the same way: 



■^ ^{diP'jif~d)+P\{f-d)dj)+^C 



Qij{r) = Q',j{T-d) - -idfP'jir~d)+P',{r-d)dj) + -dfdjp'ir-d) (24) 



and then up to quadrupole order: 

Qu = Q'.j - \ Wi + P'idj) + \didjp' + ■ ■ • (25) 

Although the individual multipole contributions do transform with a shift of 
the origin, according to (l20b (l23t (l25t and so on, the total macroscopic charge 
density: 

(?] ) = p - d,Pi + dijQ.j + --. = p'- d,P', + d.jQ'.j + ■■■ (26) 

does not, as can readily be verified. A similar conclusion can be drawn with 
respect to the macroscopic current density (j), which is given by: 

(ji) = drP, - dtdjQj, + CijkdjMk - dtP', - drdjQ'ji + etjkdjM'k (27) 

where df is the partial derivative with respect to time and where the magneti- 
zation density transforms according to: 

M,{r) = M\{r-d) - ^e,jkdjdtP'k{f-d) (28) 



and then up to quadrupole order: 

M,=M',-^-eijkdjdtP'k (29) 

Since at the macroscopic level only (t]), (j) have physical meaning, whereas 
in general p,Pj,Qij,Mi,- ■■ do depend on the choice of reference point the lat- 
ter must be considered as equivalent representations of the same physical macro- 
scopic charge and current densities. 



4 Polarizabilities 

The constitutive equations of a medium express the response fields Pi, Qij^Mi, ■ ■ ■ 
as a function of the fields Ej^Bi. This dependence can be local in space and time 
or more in general the response fields can also depend on the values of the fields 
in nearby points and/or in the past. We will use the same linear expressions as in 
Q (with slightly different notations) and we will consider only a non magnetic 
medium and again limit ourselves to terms of electric quadrupole/magnetic dipole 
order: 

Pi = P/"^ + UijEj + a.jkdjEk + G.jdtBj 

Qij = Qf+ak,jEk (30) 

Mi = MP-Gj,d,Ej 

The values of the multipole densities in the absence of any fields are indi- 
cated with a superscript (0) . In the origin dependence of the polarizabilities 
aij,aijk, Gij, ■ ■ ■ have been found based on a microscopic (and quantum mechani- 
cal) theory. Strictly speaking one should then consider first the local fields in (l30b 
and then eliminate these so that only the macroscopic fields remain. We avoid 
this complication by considering (BUI) as pure macroscopic equations where the 
polarizabilities are phenomenological parameters. From the transformation prop- 
erties of the response fields which have been found in the previous section we can 
deduce those of the polarizabilities. 

Consider first the 2nd equation which we write in extenso as Qij{r) = Q^j (r) + 
akijEj^{r). A similar equation holds for the shifted reference point Q' ijij) — 
Q' ii v) + '^' kijEk{r) . In view of (0?b we combine these 2 into: 



ij 



V /^^_/i(0)/^l:7^ n'W 



Qij{r + d) - Q',j{r) = QfMr + d) - Q']'^^?) +akijEk{r + d) ~a'kijEk{r) (31) 



Using (l24b twice this becomes: 

-\d,{P'j{r) -Pfi?)) - i(/",(F) -Pf{r))dj = akijEj,{r + d)-a'kijEk{f) 

(32) 

On the left side we use the first equation of (BUt only retaining terms up to 
quadrupole order (for the same reason the d on the right side can be dropped) and 
we then find: 

akij = a'kij - 2 (^/«'ji + cc'ikdj) (33) 

In exactly the same way we find from the second equation in (l30b and using 



Lastly from the first equation in (l30b and using (l22t we find at first: 

aijEj{r + d)+aijkdjEk{r + d)+Gijd,Bj{r + d) = a'ijEj{r) +a'ijkdjEk{f) +G'ijdtBj{r) 

(35) 

Equating the terms of dipole-order on both sides one finds: 

ccij = cc'ij (36) 

All terms of quadrupole order cancel each other if one uses (l33l . (l34b and 
also takes into account Faraday's law dtBj{r) = —ejkidj^Ei{T) and the symmetry 
of ttij. The results (l33t (l34l) and (l36b are exactly the same as those found in [i2|, 
although they have been derived in a different way. Note that in transforming 
the polarizabilities we have taken into account that with the shift of the origin 
(for calculating the moments) one must also change the (position of the) driving 
field (see e.g. (l35t . where only the first (r + d) expression is relevant due to our 
limitation to quadrupole order). It follows that if we use the origin dependencies 
of the polarizabilities as used in [Ol and combine these with the proper shift of 
the field-point then the resulting macroscopic charge and current densities do not 
depend on the choice of origin as shown by (l26b and (l27t . 

5 Reflection coefficients 

Turning now to the wave propagation in this medium it has been noted [T| that 
the (plane wave) modes are origin-independent. This was to be expected based on 



the analysis given above. In order to determine the reflection coefficients for an 
air/medium interface we need the boundary or jump conditions at this interface, 
taking into account the quadrupole effects. Choosing the (x,y)-axes into the plane 
of the interface and with the z-axis pointing into the medium these can be written 
as: 

EqAsE, = Ks = -Pz + d^Q^ + dyQ,y + djQj, (38) 

-^^ = ~Ky = -M, + dtQ^ 



A,5, = (40) 

where A, stands for the field on the air side minus the field in the medium. Tls 
is the (bound) surface charge density, K^.Ky are the components of the (bound) 
surface current density and ,^x, ^yi ^z those of the surface polarization density. 
These equations are the same as those used in except for (l38t where the contri- 
bution dvQy^ in the last term on the RHS is missing in (and the conditions differ 
from those published previously |6|). The quadrupole contributions in (l37t - (l39t 
can be understood as follows: in the bulk of the medium the quadrupole density is 
equivalent with a polarization density P^ = —djQji. At the interface this density 
becomes singular and gives rise to a surface polarization density <^, = —Qzi- The 
normal component .^^ = ~Qzz gives rise to a local voltage difference over the 
boundary and its variation along the boundary enters Faraday's law and explains 
(l37t . The in-plane components on the other hand give rise to a surface charge den- 
sity — <9x<^x — dyl^y, which explains the first 2 quadrupole contributions in (l?Sl . 
and a surface current density with components dt^x and dt^y, which explain the 
contributions in (l39b . In the bulk the polarization P^ = —djQjj is equivalent with 
a charge density —diP^ = dfdjQji and a current density d^Pf^ = —dfdjQji. The 
former has a singularity at the surface leading to a surface charge density djQj^, 
which is the last contribution in (l38t . Note that without this term the normal com- 
ponent of the accompanying current density is not balanced at the surface. 

Using Faraday's law on both sides of the interface it is easily seen that (l37t 
implies (HUb . Similarly using instead Ampere's law one can also verify that (0^ 
implies (l38t . Therefore we only need the equations (l37t and (l39b to find the re- 
flection coefficients. In general the quadrupole contributions in these equations 
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will change if another reference point is chosen, as in the bulk of the material. 
However just as in the bulk this change will be balanced exactly by the displace- 
ment of the (lower order) dipole moment and therefore if calculated correctly no 
change in reflection coefficients should be found. However whereas in the bulk 
the latter effect is accounted for automatically (and yields an origin independent 
wave equation and origin independent modes) at the free surface it must be taken 
into account explicitly. It seems indeed logical that when moving the (atomic) 
reference point, and therefore all multipoles, over a distance d, then one should 
at the macroscopic level, move the boundary surface over the same distance with 
the immediate conclusion that only d^ will be relevant. The effect of such a shift 
(alone) on the reflection coefficients is easily found to first order as: 

(J7?),, = -2Jhd,R (41) 

where "ss" stands for "surface shift" and with k^ the normal wave vector com- 
ponent in free space (we use the time dependence e-'®')- If the reflection coeffi- 
cients are origin independent then the change in the reflection coefficients due to 
the change of the quadrupole moments alone must exactly be opposite to the one 
ingB: 

(Ji?)q, = 2jk,d,R (42) 

where "qc" stands for "quadrupole change" and this should be true irrespec- 
tive of any symmetry properties of the medium. Proving this relation by direct 
analytical calculation turns out too complex so far in the most general case, this 
means for oblique incidence on an arbitrary medium ^ However we could prove 
relation (l42b if we relax one of these conditions: (i) perpendicular incidence on an 
arbitrary medium or (ii) oblique incidence on a medium with a 4-fold symmetry 
axis perpendicular to the surface. Further evidence for the general validity of (|42b 
can be gathered directly from the relevant boundary conditions (l37t and (l39b . Us- 
ing (l25t and (l29t the changes in ^^^K^^Ky due to the change of the quadrupole 
moment are given by^: 

(rf^,)q, = d,P, 

(dK,)^^ = d.dtP, (43) 

idKy)^^=d,dtPy 

Since the bulk current density has components dtPx, dfPy along the free surface 
these variations are indeed as if the surface moved outwards over d^ whereas dv > 

'"arbitrary" means in this context "without any special symmetry properties". 
^Note that due to our definition of d (R = R +d) the difference equals unaccented quantities 
minus accented quantities. 

11 



actually represents an inward movement of the surface, since the z-axis points 
into the medium. Finally we note that the lowest order term shown in (H^ is 
imaginary since it corresponds with a pure phase shift. In ||2| the variation of 
the modulus of R only was considered and therefore the relevant lowest order 
contribution to dR was missed. The variation of the modulus is indeed of a higher 
than quadrupole order and since calculations have been done up to quadrupole 
order only this is a meaningless result. 

6 Transformations 

In the frequency-domain the constitutive equations can generally be expressed as: 

(44) 



Di = £,jEj + rijBj 



where the material constants can be calculated from (l30b : 

bjj t,QOij -j- (Xij -T ] i^m\^ jmi ^imj) 

jij = jcoGij [5ij = jcoGji Vij = jIq 1 5, 



(45) 



ij 



Due to (l33t . (l34b these material constants will usually be "origin dependent" 
and therefore at first sight also material dependent properties like e.g. reflection 
coefficients. As explained in the previous section the change of reference point 
should be accompanied by a shift of the free surface boundary and when both ef- 
fects are taking into account the reflection coefficients are invariant. Raab and De 
LangeQ |l3l ||71 did not take the latter effect into account but instead, in order to 
get rid of the "origin dependence" of the material constants in (R3b . they devel- 
oped a "transformation theory", which essentially applies changes AG/y and Aa,jyt 
to the polarizabilities in such a way that Maxwell's equations remain invariant. 
The new formulation is supposed to be equivalent with the original one, but since 
both formulations lead to different physical results (in casu different reflection co- 
efficients) this cannot be true. This is due to the fact that no attention was given 
to the boundary conditions. In fact it is almost inescapable that different reflec- 
tion coefficients (in casu origin independent ones) could only have been obtained 
using this procedure by changing the boundary conditions. 

In order for the boundary conditions (l37t and ( l39b to be invariant the RHS 
terms should not change by the transformation, since the fields E and B are unique, 
thus: 

A2^, = ^K^ = AKy = (46) 
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With (l30b and since these conditions should hold for an arbitrary electric field 
it follows that: 



~AGjy + Aaj^ = (47) 

AGjx + Aajzy = 

for any/ Since in addition boundary conditions should remain invariant along 
any possible boundary (not only the (x,j) -plane considered in this example) these 
conditions should also hold after cyclic permutation x ^ y ^ z. It then readily 
follows that AGij — and Aaijk = meaning that no transformation leaves the 
boundary conditions invariant. 

7 Conclusions 

Changing the point of reference for calculating the multipole moments for the 
building blocks of a medium has two consequences. The moments themselves 
usually change with the reference point and the position in space of these elemen- 
tary multipoles changes. If both effects are taken into account then the resulting 
charge and current densities are independent of the reference point. This has been 
shown at the microscopic and at the macroscopic level. The same interpretation 
can be applied to e.g. the reflection coefficients for a semi infinite medium: if with 
a change of origin the boundary of the medium is shifted accordingly, then the re- 
flection coefficients are invariant. If the shift of the boundary is not taken into 
account then the reflection coefficient will show a variation corresponding with 
a pure phase shift. There is no need to remedy this "origin dependence" effect 
with a "transformation theory". In fact there is no transformation which leaves 
all possible boundary conditions invariant. We have performed all calculations 
including the multipoles of electric quadrupole/magnetic dipole order and for a 
non-magnetic and non-absorbing medium but we have no doubt that these results 
are still valid for more general materials. 
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